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Abstract 



We adapt Fluorescence Correlation spectroscopy (FCS) formalism to the 
studies of the dynamics of semi-flexible polymers and derive expressions 
relating FCS correlation function to the longitudinal and transverse mean 
square displacements of polymer segments. We use the derived expressions 
to measure the dynamics of actin filaments in two experimental situations: 
filaments labeled at distinct positions and homogeneously labeled filaments. 
Both approaches give consistent results and allow to measure the temporal 
dependence of the segmental mean-square displacement (MSD) over almost 
five decades in time, from ~ 40fis to ~ 2s. These noninvasive measurements 
allow for a detailed quantitative comparison of the experimental data to 
the current theories of semi-flexible polymer dynamics. Good quantitative 
agreement is found between the experimental results and theories explicitly 
accounting for the hydrodynamic interactions between polymer segments. 

Key words .-polymer dynamics, fluorescence correlation spectroscopy, F- 
Actin, semi-flexible polymers 



FCS of segmental dynamics in F-Actin 



2 



Introduction 

Living cells have remarkable mechanical properties which enable them to 
move, divide and respond to external stresses. These properties are mainly 
attributed to the dynamical characteristics of the cell cytoskeleton, which is 
a complex three-dimensional network of protein filaments mostly comprised 
of F-actin and microtubules. Both types of filaments are the polymerized 
forms of monomeric protein subunits: globular actin (G-actin) and tubulin, 
respectively. The cytoskeleton derives its strength from the elastic properties 
of these biopolymers, which, unlike synthetic polymers, are characterized by 
a high bending rigidity. 

The polymer rigidity is described by a persistence length l p above which 
thermal fluctuations can efficiently bend the polymer. Semi-flexible poly- 
mers such as F-actin and microtubules have a long persistence length: l p ~17/iin 
for F-actin (0) and several millimeters for microtubules, orders of magnitude 
larger than those of synthetic polymers. The large rigidity of these biological 
polymers enables us to experimentally address the fundamental questions of 
polymer dynamics at the sub-persistence length scales. The dynamics and 
the mechanical properties of F-actin and microtubule networks were studied 
by different opticaHechniques, such as dynamic light scattering (|2j), fluores- 
), diffusive wave spectroscopy (DWS) and microrheology 



One of the most detailed features of polymer dynamics accessible in 
experiments is the kinetics of monomer motion measured by the temporal 
dependence of monomer's MSD (r 2 (t)}. For length scales below the per- 
sistence length, the monomer displacements are anisotropic with a major 
contribution coming from the transverse modes (displacements perpendicu- 
lar to the polymer contour) and only a minor contribution from longitudinal 



moc 



to tne polymer contour) ana only a minor contribution irom longitudinal 
des. Theories HE E3, EI El, EI El, El, E3) predict that the kinetics of 



transverse motion, and thus the overall monomer's MSD, should follow the 
power law dependence of (r 2 ) oc i 3//4 . 

A feature related to the monomer MSD, the time-dependence of the lon- 
gitudinal fluctuations was obtained by measuring the end-to-end distance of 
individual actin filaments visualized by fluorescence video-microscopy (jj). 
Although the data are consistent with theoretically predicted dependence, 
the temporal resolution of video-microscopy limits the measurements to time 
scales larger than 80ms (0) and thus the overall range of (r 2 ) oc t 3 / 4 depen- 
dence spans over one order of magnitude in time only. 

A wide range of time scales was assessed using DWS and microrheology 
of micron-sized beads inserted into the F-actin mesh |8j). The (r 2 ) cx i 3//4 



FCS of segmental dynamics in F-Actin 



3 



scaling was observed from ~10/us to ~10ms. However, it is not quite clear 
how the MSD of a bead is related to the actual monomers' displacements. 
While the beads' motion depends on the properties of the mesh, the beads 
themselves may affect the dynamics of the filaments, e.g. via their large 
friction coefficient (0). In particular, the monomer kinetics measured by 
fluorescence video-microscopy (0) is two orders of magnitude faster than the 
motion of the beads measured by DWS (0). 

Here we present a new non-invasive approach to measure the monomer 
dynamics in stiff filaments: the filaments are tagged with fluorescent labels 
and the segmental dynamics of the filaments is then followed with Fluores- 
cence Correlation spectroscopy (FCS) technique. We show that FCS corre- 
lation function is directly related to the temporal dependence of monomers' 
mean-square displacement (MSD) (r 2 (t)}. The measurement of the FCS 
correlation function allows us to obtain the kinetics of F-actin monomer 
motion over a wide range of time scales, from 40/zs to 2s. Previously, 
the same method was used to study monomer dynamics in DNA polymers 
However, while the standard formulas of FCS can be applied to 
double-stranded DNA, new expressions relating FCS correlation function to 
monomer MSD have to be derived for stiff filaments such as F-actin and 
used to analyze experimental data. The main reason for that is the large 
difference in the persistence lengths of DNA and F-actin. In general, the 
dynamics of a semiflexible polymer is anisotropic: the segmental motions 
transverse to the filament are larger than longitudinal displacements. This 
anisotropy is lost at the length scales above polymer persistence length. For 
DNA, lp ~ 50nm is much smaller than typical dimensions of FCS sampling 
volume (~ 500nm) and the expressions implying isotropic dynamics can be 
used in all of the dynamic range probed by FCS. For actin filaments with 
l p ~ 17 fj,m the situation is reversed and the anisotropy of segmental motion 
has to be taken into account explicitly. 

In the next section we derive the FCS expressions for anisotropic motion 
of stiff filaments. Although the segmental motion is dominated by the trans- 
verse component, for the sake of generality, we will derive the expressions 
which take into account explicitly both transverse and longitudinal displace- 
ments. Then we use these expressions to analyze the results of two sets ex- 
periments with different labeling strategies: 1) Partially (or locally) labeled 
filaments (obtained via polymerization of nonfluorescent monomeric actin on 
fluorescent seeds), and 2) homogenously labeled filaments (obtained by poly- 
merization from a mixture of labeled and unlabeled actin monomers). The 
analysis of the experimental data with the expressions appropriate to each 
of the cases gives consistent results on temporal dependence of monomers 
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MSD. Finally, the experimental data are compared to theoretical predictions 
for the dynamics of semifiexible chains. The results are in qualitative and 
quantitative agreement with the theories taking into account hydrodynamic 
interactions between the polymer segments (fl2l . Il4h. 



Theory 

FCS technique 0, EH, (reviewed in e.g. 0, 0) is based 
on monitoring fluctuations 5I em (t) = I em — (I em ) in fluorescence emission 
I em {t) as fluorescence species diffuse in a spatially restricted excitation field, 
formed typically with the help of confocal optical scheme |2?h. The auto- 
correlation function G(t) = (5I em (0)5I em (t)) of emission fluctuations reflects 
the kinetics of motion of fluorescent sources. 

In this section we adapt the general formalism of FCS to the case of 
dynamics of linear stiff polymers and derive expressions relating FCS cor- 
relation function to the temporal dependence of MSD (r 2 (t)} of polymer 
segments. 

The instantaneous detected emission, average emission and the corre- 
lation function of fluorescence fluctuations are found through the spatial 
distribution c(r, t) of fluorescent labels and excitation-detection profile 1(f) 

hm(t)=Q f dfl(f)c(f,t), 

(Iem) = QcJdrI(f), (1) 

G{t) = Q 2 J drdr'I{r)I(f')(5c(r, 0)Sc(f', t)), (2) 

where Q is specific brightness of a fluorescent molecule dependent on fluo- 
rophore properties and the efficiency of detection optics, c = (c(r, t)) is the 
average concentration of fluorophores, and 5c(r, t) = c(r, t) — c. Following 
13) Eq. 12 can be rewritten as: 

G(t) = J dq\I(q)\ 2 (5c*(q,0)Sc(q,t)), (3) 

where "tilde" denotes spatial Fourier transform of the corresponding quan- 
tities, such as e.g. I(q) = (2ir)~ 3 / 2 J df I[f)e~ tqr , and V is the total volume 
of the sample. 

We assume that fluorescent molecules are not moving independently, but 
are attached to relatively large objects which themselves move independently 
of each other. We will further assume that the objects are statistically 
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equivalent and that the distribution of fluorophores within each object with 
respect to its center of mass is described by 3>(r, s(t)), where s(t) denotes 
the set of internal degrees of freedom describing the current conformation 



the center-of-mass position and internal conformation of object j at time t. 
Finally, substituting these formulas in Eq. |3]we have: 



where n is the average concentration of the objects (number of objects per 
unit volume), Ar(i) = fj(t) — Vj(0) is the displacement of an object, and the 
index j was omitted everywhere in Eq. 0] due to the statistical equivalence 
of the objects. 

In its general form Eq.|3]can be applied to any objects which have some 
internal structure and internal dynamics. 

We assume now that the objects are uniformly labeled segments of semi- 
flexible polymers (one labeled segment per polymer, which, in general, are 
longer than their labeled parts and contain unlabeled parts) . For sufficiently 
stiff polymers, we can neglect the dynamics within the labeled segments and 
consider them to be straight. In this case the set of internal degrees of free- 
dom s reduces to a unit vector s defining the orientation of the segment. 
Within the same approximation we can assume that any given labeled seg- 
ment moves without change in its orientation s)(0) = Sj(t). We will discuss 
in detail the validity of our assumptions at the end of this section. Here we 
just note that although these assumptions may look prohibitively restrictive, 
for sufficiently stiff filaments they are valid in a wide range of segmental dis- 
placements. 

For a thin straight segment of length L uniformly labeled with linear 
density a of fluorophores: 



where a is an angle between q and s. Furthermore, we split the segmental 
displacement into components parallel and perpendicular to the segment 
Ar = Arii + Ar_j_, we assume the two components to be independent of each 
other and to be Gaussian random variables. Then, for fixed q and s the 





(4) 
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ensemble average of e is given by: 

q 2 cos 2 a , . 9 . q 2 sin 2 a , . 9 . \ 

-^—2 — < A i> - ^—4 — < Ar ' ) ) ( 6 ) 

The difference in the numeric pref actors in the two terms in the Eq. cisterns 
from the fact that Arn is defined on a line (parallel to s), while Afj_ is 
defined in a plane (perpendicular to s). 

We assume, as it is usually done for confocal setups, the excitation- 
detection profile I {f) to be 3D Gaussian axisymmetric with respect to optical 
axis Z: 

/(r) = /oexp(-H(£W)-g) (7) 

I(q) = exp (-^V sin 2 B - fq 2 cos 2 0) , (8) 

where w xy and w z define the width of the profile in the XY plane and in Z 
direction respectively, and 6 is an angle between q and Z-axis. 

Substituting © , © and (JBJ) into Eq. QJ averaging over a for a given q 
and integrating over q, we have: 

G(t) = ^I 2 wl y u 2 Q 2 a 2 n dk j\u (^) exp {-^f{h\,h 

(9) 

where u = cos 9, p = cos a, reduced units k = qw xy , X = L/w xy , /i 2 = 
(Ar^/w^y, h 2 , = (Ar 2 )/w xy and to = w z /w xy were introduced, and / de- 
notes the following expression: 

f(h 2 ± ,hj,p,u) = h 2 ± -p 2 (h 2 ± -2hj) + u 2 (uj 2 -1) + 1 (10) 

FCS correlation function is usually normalized by the square of the av- 
erage emission G\(t) = (5I em (0)SI em (t)} / {hm} 2 • With this normalization 
the amplitude of the correlation function at short time scales is the in- 
verse of the average number N of molecules in the detection volume (given 
by ir 3 / 2 w xy w z ): G\(t — > 0) = 1/N. Here we prefer another normaliza- 
tion: G 2 (t) = (5I ern (0)SI em (t)) I (I em ) = (hm)Gi(t). A correlation func- 
tion defined this way is independent of the concentration of the moving 
species (as long as there are no interactions between the objects) and its 
amplitude at short time scales gives the fluorescence per moving object: 
G%[t — ¥ 0) = (I em )/N. This is an interesting quantity in the context of 
labeled segments comparable or larger than w xy : in this case only part of 



(e-^ Ar(t) ) k> -=exp 
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the labeled segment can "fit" into the detection volume and contribute to 
the correlation function. The length of this part can be estimated from 
G 2 (t-0). 

In order to calculate G2(t) we substitute c = nuL and (J7J) into (pQ) to 
find 

(I em ) = (<K/2f/ 2 I w 2 xy w z QnaL. (11) 

Performing integration over k in Q and making use of (jllj) and of definition 
of G2 we obtain: 

IoQ w y 1 , /-i l-exp(-AV//(^,^,p,n)) 

G2U) = — ~r~ / an / an , . (12) 

2{> 2V2 xy XJo Jo P i?JHhl,hlp,u) 



Given the knowledge of the experimental geometry (i.e. parameters w xy , w z 
and L), Eq. El can be used to numerically calculate the relation between 
the FCS correlation function and temporal behavior of segmental MSDs, i.e. 
h\(t) and h?,(t). 

In two important limiting cases the explicit expressions can be derived 
for G2(t): 1) for very short labeled segments, i.e. L — > while La = const, 
and 2) for very long labeled segments, i.e. L — » 00 while a = const. 

In the case of short segments, we obtain: 

nM IoQ *Lu , 

G2\t) = — p — - = In ■ 



2v / 2^(l + / l 2 L )(o;2_ 1) ( /l ^_2/ l 2) " y J(l + 2hp(u 2 + h 2 ± 

(13) 

We note that for isotropic motion, i.e. for h\ = 2h 2 = (Ar 2 ) / (3w 2 y ) , 
the Eq. ED reduces to the more standard FCS expression for the random 
motion of point-like objects: 

For infinitely long segments Eq. 1121 gives: 



(h(t) = 91° ™*«»^ arc tan Jp£. (15) 



Note that ^(t) in this case is independent of /i| since longitudinal mo- 
tion of infinitely long labeled segments does not lead to any fluctuations in 
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fluorescence. 

Although unrelated to our experiments, still an interesting particular 
case of application of Eq. ^] is that of spherically symmetric detection vol- 
ume w xy = w z = w. In this case (^(t) = QIowa(l + h 2 ^)^ 1 is similar 
to the correlation function produced by a planar motion of point-like objects. 

We return now to Eq. ^] in order to find the dependence of the flu- 
orescence per moving object G%(t — ► 0) on the length of the labeled seg- 
ment. We make use of the fact that for small segments L <C w xy : (^(t — > 
0) = QIq(tL/(2v2), and define an apparent length L app and respectively 
A a pp = L app /w xy such that for the labeled segment of any length G^it — ► 
0) = QI aL app /(2V2). 

For t = 0, i.e. /i| = fij_ = it is possible to perform integration over u 
in Eq. El an d arrive to the expression relating the apparent length of the 
labeled segment to its actual length and to the parameters of the detection 
volume: 



A 



\n(co + \Juj 2 - 1) + f — £= ( Ay^erf ( - ] + ve~ 
Ji vVv 2 - 1 V \ v ) 

(16) 

where erf(x) = 2ir 1 / 2 e 1 dt 

An example of A app (A) dependence is shown in Fig. ^for uj = 5 corre- 
sponding to our experimental geometry. As expected for A < 1, all of the 
segment can fit into the detection volume and X app ~ A). For A > 1 the 
apparent length saturates at the value which can be found from 



Aapp(A — > oo) — UaI — 2 7 arctan Vto 2 - 1. (17) 

Finally, we discuss the validity of assumptions leading to Eq. 1121 The 
main assumption we made was to neglect the internal dynamics within the 
labeled segment. This assumption in fact just puts a lower limit on the 
accessible range of studied segmental displacements: the derived equations 
are valid as long as the center-of-mass motion Af(t) of the segment is larger 
than the characteristic motions within the segment. 

The characteristic internal motions (Arf nt ) within the labeled segment 



can be estimated to be of the order of {Ar 2 nt ) = (2/45)L 3 // p (14). For 
segmental displacements larger than that, motions within the labeled part 
can be neglected and the labeled segment can be considered essentially rigid 
and moving as a single unit. 

This condition seemingly prohibits studies with long labeled segments, 
which have considerable motions within them. However, only a small part 
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of the filament (of length ~ L app < L) can cross the sampling volume and 
contribute to the fluctuations in fluorescence at any given moment. Thus 
the lower limit for the range of accessible segmental motions can be further 
relaxed to (2/45)L^ pp // p . E.g. even for very long homogeneously labeled 
actin filaments (L S> w xy , L app 2.5w xy ) the segmental dynamics can be 
studied in the range (Ar 2 ) > 4 • 10 -4 //m 2 . 

The expressions derived in this section (such as Eqs. El DID Ej) can 
be directly used to measure segmental displacements from FCS correlation 
functions: only the parameters defining the experimental geometry (w xy , w z 
and in the case of Eq. ^] L in addition) need to be calibrated. All other 
parameters affect G%(t — > 0), the value of which can be determined from the 
plateau level of the experimentally measured correlation function at short 
time-scales. Finally, for the case of internal dynamics of a semi- flexible 
chain, relevant to most of the studies on actin dynamics as well as to our 
studies presented here, the longitudinal displacement h?, can be neglected 
in comparison to transverse motion h\. Then, with known geometrical pa- 
rameters and measured G^t — > 0), Eq. ^Jand its limiting cases Eas. Upland 
1151 give a one-to-one relation between segmental MSD and FCS correlation 
function. Some examples of such dependence are given in Fig. [21 Using 
these dependencies, the experimentally measured correlation function G^t) 
can be converted into the temporal dependence of segmental MSD (Ar^(t)} 
{=wl y h\{t)). 

Materials and Methods 
Actin preparation 

Unlabeled actin is purified from chicken skeletal muscle acetone powder and 
stored in G-buffer |29l). Two sets of samples are prepared: 1) filaments 
labeled at defined positions, and 2) homogenously labeled filaments. 

In order to label the filaments at defined positions (sample 1), we utilize 
fluorescently labeled actin filament seeds as templates for additional poly- 
merization of unlabeled actin monomers. To prepare fluorescent seeds we 
first polymerize 1/xM of fluorescent G-actin (Actin Alexa 568, Molecular 
Probes, Eugene, OR, or Rhodamine-Actin, Cytoskeleton, Denver, CO) in 
the presence of phalloidin (Molecular Probes) to stabilize the filaments (1:1 
actin to phalloidin molar ratio). The labeled filaments are then broken by a 
brief sonication and vigorous pipetting into short fragments (average length 
of ~ 170nm, estimated by FCS see Results and Discussion). These frag- 
ments are then used as seeds for further polymerization of unlabeled G-actin 
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(9/iM). Polymerization proceeds for 10 min at room temperature. We note 
that concurrent with polymerization, there is an ongoing annealing process 
which both increases polymer length and creates multiply labeled fil- 
aments. Independent fluorescence microscopy observation confirms that at 
the beginning of the experiment about ten percents of filaments are labeled 
at two distinct positions (typically separated by more than 1/im), while the 
rest are single-labeled. Since the distance between the labeled portions of the 
double-labeled filaments is much larger than confocal radius (0.21 /xm), the 
"cross-talk" between the labels can be neglected and the formalism derived 
in the preceding section can be applied. 

The homogenously labeled filaments (sample 2) are prepared in a simi- 
lar manner. The difference is that the seeds are prepared from a mixture of 
labeled and unlabeled actin monomers (1:9 molar ratio) and further poly- 
merization proceeds with the same mixture. 

This procedure results in actin filaments of several microns in length: 
~ on average at the beginning of the experiment and growing in the 
course of experiment to ~ 8/im due to annealing, as verified by fluorescence 
microscopy. 

For the experiments, the solution is diluted tenfold to a final actin con- 
centration of 1/xM. Typically, one microliter of solution is sealed between 
two glass coverslips separated by a 250/um spacer. To prevent protein ad- 
sorption, the glass coverslips are coated with an inert polymer (Polyethylene 
Glycol) according to the protocol of (31). Most of the measurements were 
carried out at a distance of 40/im from the surface. This distance was chosen 
on purpose to be larger than F-actin length in order to minimize the effect of 
surface proximity on filament dynamics. Control experiments performed at 
a distance of lOO^m from the surface give results identical to those presented 
here. 

The experiments are started immediately after dilution and are con- 
ducted within 30 minutes. The main reason to limit the duration of experi- 
ment is to minimize the effect of actin filament fragmentation, which leads 
to appearance of short filaments and associated noise in the FCS correlation 
function. After ~1 hour of FCS measurement we start to see the changes in 
the MSD due to the fragmentation: the MSD vs. time curve shifts to larger 
displacements. To be on the safe side, we limit all of our measurements to 
the first 30 min after polymerization. 
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Experimental setup 

The optical setup is home-built based on the Nikon Eclipse TE300 inverted 
microscope (Nikon Corporation, Tokyo, Japan). The confocal excitation is 
provided by 514nm line (~ 2.5^W power before microscope objective) of 
an Ar-ion laser (Advantage 163D, Spectra-Physics, Mountain View, CA) 
deflected by Q525 dicroic beamsplitter (Chroma Technology, Rockingham, 
VT) into a high-power objective lens (UPLAPO 60X1. 2 W, Olympus Europe, 
Hamburg, Germany). The collected emission passes through the beamsplit- 
ter, then a bandpass filter HQ565/80 (Chroma Technology, Rockingham, 
VT) and a pinhole of 50fim in diameter. The emission is detected by a 
photon counting avalanche photodiode (SPCM-AQR-14 PerkinElmer Op- 
toelectronics, Vaudreuil, Quebec, Canada) whose output is fed into digital 
correlator Flex2k-12Dx2 (Correlator.com, Bridgewater, NJ). The correlator 
is capable of working in two modes, either as traditional correlator carrying 
out the correlation analysis of emission online, or as photon history recorder, 
storing the time arrivals of every photon on computer hard drive. For the 
presented experiments, we make use of the photon history recorder mode, 
while analyzing the recorded photon traces offline with software correlator 
as described in Data Analysis. The correlator program was written as C- 
module running under MATLAB environment (MathWorks, Natick, MA). 

The parameters w xy ~ 0.21/im, w z ~ l.lfim of the confocal volume are 
calibrated before and after each experiment by measuring the diffusion of 
free Rh6G fluorophores (0). 

Data Analysis 

Sample 1: The measurement of the photon emission count rate from lo- 
cally labeled F-actin reveals that photons arrive in intense bursts of ~ 10 5 
counts/sec lasting ~ 0.1 to Is, separated by intervals of low count rate of 
~ 10 3 counts/sec (Fig. |SJ). The bursts are caused by the passage of the 
labeled F-actin through the confocal volume. The fluorescence in between 
the bursts originates from residual free fluorophores diffusing in the sample. 
The motion of the free fluorophores results in a correlated background noise 
I b (t) which adds up with the labeled actin signal I(t) to a total emission 
hot = I+h- Thus the overall correlation function Gtot{t) = {8Itot(0)8It o t(t)) 
is given by: 

G tot (t) = G(t) + G b (t), (18) 

where G(t) = (8 1 (0)5 1 (t)) and G b (t) = (8I b (0)8I b (t)) are the correlation 
functions of emission from labeled segments and from free fluorophores re- 
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spectively. 

Although the overall contribution of free fluorophores Gf, to the total 
correlation function is small, their fast motion is responsible for most of the 
decay of the correlation function at the time scales below 1ms, where the 
motion of the filament segments is negligible. This could limit the analysis of 
monomers' MSD to the time scales above 1ms. However, as shown below, it 
is possible to separate the contributions of labeled actin and free fluorophores 
within the same experiment, and thus, extend the range of measurements 
to time scales as low as ~ 40/is. 

In order to separate free fluorophore noise from the signal, we record 
the complete photon trace, i.e. times between arrivals of consecutive pho- 
tons (with temporal resolution of 16.7ns). The photons are then binned 
into 100ms intervals {I(t n )} (Fig. |2J)- The stretches of time with no bursts 
are determined and the background correlation function G&(i) is calculated 
on these stretches using the original photon traces. The total correla- 
tion function Gtot(t) is computed using the complete photon trace. Fi- 
nally, the contribution of labeled actin G{t) is obtained using Eq. ()18j) . 
The intensity-normalized correlation function of actin segments is found by 
G 2 {t) = G{t)/{I tot -I h ). 

We find this procedure more robust than the correlation analysis of 
bursts intervals. First, there is some background contribution within the 
bursts as well. Second, unlike bursts, the intervals of pure background can 
be determined unambiguously: any intervals suspect to contain a burst can 
be deselected. 

Practically, the background is deduced by calculating the median inten- 
sity over all bins I me d = median({I (t n )}) , and selecting all of the intervals 
which deviate from the median by less than a = median({(I(t n ) — Imed) 2 }) 
for analysis of background noise. An example of this procedure is presented 
in Fig. 01 The resulting total, background and signal correlation functions 
Gtot{t), Gb(t) and G(t) are shown in Fig. ^ The background correlation 
function Gb{t) is indeed well described by 3D diffusion model of the free 
fluorophore 1)271 ) with a characteristic decay time of 70/iS. We note that 
the amplitude of G(t) is much larger than that of Gb(t) mainly due to the 
multiple labeling of the filament segments: thus the motion of the labeled 
segments leads to much larger fluctuations in emission than the motion of 
single fluorophore molecules. 

For calculation of monomers' MSD the amplitude of the correlation func- 
tion Go = G(t — > 0) is estimated from the level of the correlation function 
in 3 to 20 //s range, which is above the characteristic time scales of the 
fluorophore triplet state kinetics and below the characteristic time scales of 
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segmental dynamics. 

Sample 2: We use the same approach to analyze data from homoge- 
neously labeled F-actin. However, since the passage of the labeled segments 
through the sampling volume is more frequent in this compared to 

Sample 1, we make use of shorter binning intervals of 30ms and we pick 
the intervals with the average count rate not exceeding I me d + 0.25<r for 
the analysis of background noise. These conditions give a noise correlation 
function with characteristic decay time below 100/US. The above parame- 
ters were found to be optimal between less restrictive conditions which lead 
to a notable contribution of the signal in the estimated noise correlation 
function (characterized by decay times exceeding 1ms), and more restrictive 
parameters which clearly underestimate the noise level. 

Results and Discussion 

We present the correlation functions obtained from both types of samples in 
Fig. 21 To facilitate the comparison of temporal kinetics, the amplitude of 
the correlation function of homogeneously labeled actin was adjusted to the 
level of the correlation function of the partially labeled sample. The func- 
tions look similar but are notably shifted in time: the correlation function 
of the homogeneously labeled sample decays slower by a factor of ~ 1.4 than 
that of the locally labeled F-actin. 

We can estimate the length of the labeled parts of Sample 1 by analyzing 
the amplitude of its intensity- normalized correlation function (^(i — > 0). 
Ct2 amplitude is 62 ± 8 times higher than the corresponding amplitude of 
the correlation function of G-actin monomers obtained in similar conditions 
(data not shown). Since 14 actin monomers form a filament of 37nm, this 
gives the apparent length of the labeled segment of L app as 160 ± 20nm. 
Converting L app into real segment length (Eq. ED or Fig. ^) we obtain L ~ 
170 ± 30nm. 

The apparent length of homogeneously labeled F-actin is ~ 470nm ~ 
2.3w xy , as estimated from the amplitudes of the corresponding correlation 
functions. This value is in a good agreement with our expectations for long 
homogenously labeled filaments (Eq. ^Jand Fig. 

As discussed above, in order to extract (Ar^(t)) dependence from the 
correlation functions we neglect longitudinal motion of the segments. Al- 
though the general expression Eq. ED can be used to analyze both sets of 
data, the numerical calculation shows that for L < w xy (i.e. A < 1) Eq. El 
leads to essentially the same dependence of the correlation function G2 on 
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MSD as Eg. 1131 derived for point sources (Fig. [2J • Similarly, for homoge- 
neously labeled filaments with length larger than 3/im (A > 15), the Eg. 1151 
for infinitely long labeled segment can be used (compare curves for A > 15 
and A — > oo in Fig. [2J) . Thus we make the use of explicit expressions Eq. 
and Eq. El to analyze data on partially and homogeneously labeled poly- 
mers respectively. Practically, in both cases we tabulate G^/i^) / ^2(0) for 
a wide range of h\ as exemplified in Fig. [21 (curves for A = and A — > 00). 
We normalize the experimentally obtained correlation functions G 2 {t) by 
their plateau values. We use tabulated G 2 (h 2 L )/G 2 (0) depend encies to find 
h\ values corresponding to each measured data point G 2 {t). This gives the 
h\{t) dependence (and respectively (Ar^t)) = w^. y hj_(t)). 

The extracted temporal dependencies of transverse monomer motion are 
presented in Fig. El Despite the difference in the temporal behavior of 
the original correlation functions for partially and homogeneously labeled 
filaments, the application of the appropriate expressions for each case leads 
to consistent data on monomers dynamics for both sets of measurements. 
The data span a wide range of time scales: from ~ 40/is to ~ 2s. 

In Fig. |5] we compare the experimental data to the predictions of two 
types of semi-flexible polymer dynamics theories, which do (flih and do not 
account for the hydrodynamic interactions (fl^h. The measurements agree 
very well with the hydrodynamic theory over about four decades in time. 
The deviation of the data at long time scales is probably related to the 
diffusion of the filament as a whole. 

To conclude, we adapt FCS formalism for the studies of the internal dy- 
namics of semi-flexible polymers. We make use of a developed formalism 
to obtain non-invasive measurements on the kinetics of segmental motion in 
actin filaments. Two labeling strategies, local labeling and homogeneous la- 
beling, lead to consistent results. The transverse segmental MSD (Ar^_(t)) 
is probed over a wide range of timescales, from ~ 40/is to ~ 2s. Almost 
over the whole range the data points follow closely the prediction of hy- 
drodynamic theories |l2l . Hih . Thus, non-invasive measurements of (r 2 (t)) 
carried out with FCS allow us to test hydrodynamic theories directly over 
a wide temporal range. We note, finally, that although in the presented 
measurements the transverse motion is the dominant mode of motion, in 
other cases, such as F-actin in the presence of molecular motors (myosins) 
the longitudinal motion may be significant. Then the two suggested labeling 
strategies, partial and homogeneous labeling, being very different in their 
sensitivity to longitudinal motion, can be used in conjunction to separate 
transverse and longitudinal components of segmental displacement. 

We are indebted to R. Granek, E. Frey and K. Kroy for fruitful discus- 
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reading of the manuscript. This work has been supported by the Israel 
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Figure Legends 

Figure ^ 

Dependence of the apparent length of the labeled segment on its real length 
as given by Eq.EDf° r w = 5 (close to the aspect ratio of the sampling volume 
in our setup). The lengths are given in the units of confocal volume radius 
A = L/w xy , X app = L app /w X y . The apparent length of very long A > 1 
segments approaches « 2.5w xy . Inset: Ratio of the apparent length to the 
real length of the labeled segment. Labeled segment lengths of L < w xy 
(A < 1) results in L app L. 

Figure 

Calculated dependencies of FCS correlation function on the transverse mean- 
square displacement of semi-flexible polymer segment. FCS correlation 
functions are normalized by their zero-time values (to have unit ampli- 
tude) and transverse MSD is given in the units of confocal volume radius 
h\ = (Ar'j) /w xy . The leftmost curve is given by standard FCS expression 
Eq. E]f° r isotropic motion. The other curves are calculated from Eas. 1121 
El El for different values of A, left to right: 0, 1, 2, 5, 15, oo. Inset: Same 
curves in semilog scale allowing to assess wider range of h\. 

Figure [HI 

Photon count trace of partially labeled actin filaments (Sample 1). The pho- 
ton trace is collected with 16.7ns resolution and split into 100ms bins. The 
data points represent photon counts per bin: (a) Full trace over the duration 
of one measurement, (b) Zoom into first 30s of measurement. Bursts (dotted 
line) due to the passage of labeled parts are separated by the intervals of 
background noise (solid line). 

Figure |3J 

Total Gtot(t) (dotted line), background Gi(t) (thin solid line) and signal 
G(t) (thick solid line) correlation functions of partially labeled filaments 
(Sample 1). Gtot(t) is obtained by analyzing the complete photon trace, 
Gb(t) is caused by the diffusion of free fluorophore and is calculated from 
the photon trace in between the bursts in Fig. El G{t) = Gtot(t) — Gb(t) is the 
correlation function resulting from the motion of F-actin labeled segments 
after background substr action. Dashed line is the correlation function of 
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homogeneously labeled actin (Sample 2) after subtraction of background 
noise. The correlation function of Sample 2 was normalized to have the 
same amplitude as that of Sample 1 in order to facilitate the comparison of 
temporal behavior. 

Figure \5[ 

The kinetics of random motion (Ar^(t)) of actin filaments' segments. Ex- 
perimental measurements on locally labeled (thick dashed line) and on ho- 
mogeneously labeled (thick solid line) are compared to the theoretical pre- 
dictions by the hydrodynamic theory of Granek (14) (thin solid line) and 
non-hydrodynamic theory by Harnau et al (flit ) (thin dashed line). The pa- 
rameters used for the calculation are l p = 17/um, filament diameter of 7nm, 
solvent viscosity ImPa ■ s, and filament length of 6/im. 
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